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, Abstract. We present calculations of Post-Newtonian parameters for Brans-Dicke 

CN| ' tensor-scalar gravity in an arbitrary number of compact extra dimensions in both 

the Jordan and Einstein conformal frames. We find that the parameter 7, which 
O ! measures the amount of spacetime curvature per unit mass, becomes a function of w, 

the coefficient of the scalar kinetic term in the Brans-Dicke Lagrangian. Experiment 
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^ . has placed strong constraints on 7 which require that uj become negative in the Jordan 

frame for any number of extra dimensions, highlighting that this formulation is not 
physical. We also confirm the well-known result that a compact extra dimension can 
be equivalently viewed as a massless scalar 'dilaton.' In the Einstein frame, we find 
' that the behavior of 7 as constrained by experiment replicates that which is predicted 

CO . by string theory. 
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' 1. Introduction 

The addition of a scalar field is one of the simplest possible modifications to general 
relativity. Tensor-scalar models have been studied as far back as 1948 (Thiry 1948, 
Jordan 1959). Theories with a massless scalar were introduced by Brans and Dicke 
(1961) to satisfy considerations related to Mach's principle. However, Brans-Dicke 
gravity continues to be of interest today as string theory also predicts a massless scalar 
component of gravity {e.g. Lidsey et al 2000). 

Another prediction of string theory is compact extra dimensions. Again, such a 
notion has been considered for many years in the context of unification, perhaps most 
notably in the classic papers of Kaluza (1921) and Klein (1926). 

It is therefore of considerable interest to study the behavior of Brans-Dicke type 
theories in an arbitrary number of compact extra dimensions. The Parametrized Post- 
Newtonian (PPN) formalism (cf. Will 1993) lends itself to the extraction of observable 
parameters in a theory-independent manner. Xu and Ma (2007) calculated the PPN 
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parameters for general relativity with a single compact extra dimension and found that, 
although theoretically viable, it is ruled out by experiment. Here we consider Brans- 
Dicke gravity in an arbitrary number of extra dimensions. 

We complete the calculations in the Jordan frame as the theory was originally 
formulated. Nevertheless, several arguments have been made demonstrating that the 
Jordan frame formulation is in fact unphysical (Faraoni et al 1999; Magnano and 
Sokolowski 1994; Cho 1992). Our results confirm this conclusion. In the following 
discussion, we then also transform our results to the physical Einstein frame, and 
determine the constraints which have been placed on the PPN parameter values by 
Solar System experiments. 

The organization of this paper is as follows. In §2 we discuss the overall approach 
to the calculation. In §3 we present the actual PPN calculations. Finally, in §4 the 
implications of the PPN coefficients are discussed. 

In this paper we follow the signature convention ( — h + + . . .). In §3 for the sake 
of brevity, the sign indicates equality to the appropriate order. 



2. Overview of the Calculation 



The Brans-Dicke model postulates a dynamical massless scalar field which couples 
directly to the Ricci scalar and also includes a kinetic term whose weight is set by the 
single parameter cj, yielding a Lagrangian density 

^^R(t>- {u/(t>)d^(t>d>'(t> + =5f^atter. (l) 

The field equations are obtained by varying (1) with respect to and the metric g^i,, 
and read 

Stt uj 

+l(V.a^0-5^.V,V''0), (2) 

^v,v''0-^(a0f = -i?, (3) 

where it! = g^uR^'^ is the Ricci scalar, (90)^ is shorthand for 9^09'^0, and is the 
covariant derivative. 

We will consider Brans-Dicke theory in d + 1 dimensions. Since we know that our 
universe contains three large spatial dimensions, — 3 of the dimensions considered here 
shall be compact. Making use of the fact that in + 1 dimensions g^j^uQ^'^ = + 1, we 
may contract the field equations in order to write them in the more useful forms 

where again T = g^^T'^''. 
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Lastly, we will also assume a perfect fluid matter source such that 

T^^- = (p + pn + p)u^u'' + pg^", (6) 

where p is rest mass energy, H is internal {e.g. thermal) energy, and p is isotropic 
pressure. 

We wish to solve the field equations for the metric to post-Newtonian order; that 
is, we want /i^,^ such that 

goQ = -1 + /loo ~ 0(4) 

Qqi = ~ 0(3) 

^,, = l + /i,,~0(2) (7) 
in terms of the d dimensional potentials defined as 
V^C/ = -47rp, VVj = -^-KpVi, 
V^$i = -A'Kpv'^, V^$2 = -47rpt/, 
V2$3 = -Anpll, V^$4 = -47rp, 

where indicates the spatial Laplacian. 

In accordance with standard PPN bookkeeping, we count velocity v <^ 1 ~ 0(1). 
Virial relations then give v"^ ^ U p/ p li 0(2). Finally, = v - V ^ 0(1) where 
V is the spatial gradient operator. 

Having obtained the metric to sufficient order, we will compare it to the standard 
expansion of the metric in terms of the PPN parameters, 

^00 = -1 + 2C/ - - 2e$i^ + (27 + 2 + + Ci - 20^1 
+2(37 - 2/3 + 1 + C2 + 0*2 + 2(1 + (3)^3 
+2(37 + 3C4-20$4- (0-20^ 

9m = -|(47 + 3 + ai - a2 + Ci - 20^i - |(1 + «2 - Ci + '^OWi 

= (l + 27[/)%. (8) 

This expansion of the metric has been specifically designed to endow these parameters 
with physical significance. The parameter 7 is seen to specify how much space curvature 
is produced per unit mass, while (3 indicates the strength of nonlinearity in the 
superposition law. Both these parameters are equal to one in general relativity, with all 
remaining parameters equal to zero. The a parameters measure preferred frame effects; 
C measures violation of conservation of energy-momentum; and ^ measures preferred 
location effects. More details can be found in a reference such as Will (1993). 

3. PPN Coefficients 
3.1. C to 0(2) 

To begin, we must solve the field equation for 0. We may expand the scalar field around 
a background value 4>o as = 0o + C + 0(4). This allows us to write 

= (00 + C)-' = 0o"'(l - C/00 + 0(4)). (9) 
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The 0(2) term will only be necessary when we are working to 0(4) overall. The only 
0(2) term in the energy-momentum tensor is = p while Qqq — —1 + 0(2) so that 
T = —p. The field equation is now 

Since do ~ 0(1), The D'Alembertian may be reduced to the usual spatial Laplacian. 
The equation may now be solved in terms of the Newtonian potential U as 

(11) 



where 



2 d-l + {d-2)u; 

00 d+{d-l)u; ■ ^ ^ 



3.2. hoo to 0(2) 

To 0(2), (4) becomes 

(13) 

which is solved by 

hoo = '2GU. (14) 
Thus we recognize G as Newton's constant which is now a function of the parameter cu. 

3.3. hij to 0(2) 

Gauge freedom in general relativity allows us to make the transformation h^i, — > 
+ d^^i, + dn^^n for an arbitrary vector field Hence by choosing to satisfy 
dfxd'^ii — —d'^h^i + + di(/(f)o we may enforce the gauge conditions 

dX - Id.h^, = i-d.C (15) 
This allows us to cast Rij in the form 

i^.^_lVX + i^.^.C- (16) 

As Tij is at least 0(4), it may be neglected, and the right hand side of the field equation 
(4) becomes 

Setting these equal, we obtain 

From here we may read off the PPN coefficient 

^= ^-±^ (19) 

^ d-l + {d-2)u;' ^ ^ 



= -2G^.:r4T7x4vT + '2Gdod,x, (23) 
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3.4. hoj to 0(3) 

Imposing one additional gauge condition 

d,K - \d^K = l^oK + (20) 
allows Roj to be written 

R^j = ^dodX + kdoid.h'; - \d^h^^) - IV'hoj + ^^dodjC. (21) 

The trace term on the right side of the field equation is now 0(5) and may be neglected 
leaving 

^o, = -^pt', + ^MC. (22) 

These equations are satisfied if 

d+{d- 1)00 
'd-l + {d-2)uj 

where Vj is the vector potential defined above and x is the "superpotential" defined by 
V^X = — 2t/. This may be rewritten using the potential Wj from Will (1992) as 

3.5. hoo to 0(4) 

We now return and extend the evaluation of hoo to 0(4). Applying our gauge conditions 
and previous results to the full expansion of the Ricci tensor yields 

D _ iv72,(4) G ^ 2ci-3 + (2(i-4)a; , , 

i?oo - -,V hoo + d-l + {d-2)J'^ - d-l + {d-2)u ^ ^^^^ 

+i—r^^^r^G'uv'u. (25) 

d-l + {d-2)ou ^ ^ 

Using the normalization condition u^u^g^^, — —1, we can solve for 

^0^0 ^ "^^u^n^grj ^^^2GU + v^ + 0(4). (26) 
-500 

Substituting into (6) gives, to the necessary order, 

Too = T°°5o^ = p(H- n - 2GU + v""). (27) 

Meanwhile, 

T = T^'^g^, = -p(l + n - dp/p). (28) 

Recall that we must now also include the factor 1 — C/0o = l — GU /{d — l-\-{d — 2)ijj) 
from (9). We also must evaluate the last term in (4) carefully. 

Vo^oC = {dl - T^MC = {dl + Gd'Udi)C (29) 
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Rearranging, we find that 



/ 00 



d - 1 + ((i-2)o;p. 

2G''UV'U, (30) 



d - 1 + (ci-2)a; 
which is solved by 

hoo = 2GU - 2G'U' + 2( | ~ 



d - 1 + (o? - 2)u; 

+2cif— — i^^^— )g$4, (31) 



.d-l + {d-2)u. 
where we make use of the fact that 

V^{U^) = 2UV'^U + 2{yUf. (32) 

3.6. Results 

We must now dimensionally reduce our results before comparing them with the standard 
four dimensional form of the PPN metric. This may be accomplished for any source 
field by integrating over the compact extra dimensions. For example, in n = d — 3 extra 
dimensions, 

= j d"a: (33) 

where g' is the determinant of the n x n extra dimension sector of the metric. In the 
post-Newtonian expansion, this part of the metric will be diagonal so 

X 1/2 

^/W\-[I[9ii) =(1 + 27C/)"/^ (34) 

Using the Green's function G, the d+1 dimensional potential may be written in terms 
of the four dimensional potential plus an 0(4) correction as follows, 

U{x) = J G(|f-f'|)p(f) d^x' d"f 

' G(|f - f |)p(^)(f )(1 - n^U)d^x' = C/(^)(f) - n7$2(f). (35) 



In fact this is the only correction we need consider. Since the integrands of other 
potentials are already 0(4), their corrections will be 0(6) and may be neglected. 

Effecting the change of variables in (31), we may now compare our results with the 
standard PPN metric (8), to find 

^= d-l + id-2p ^ ^ = = 
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with all others zero. 
4. Discussion 

We first note that there is degeneracy within 7 between d and uj. For example, 7 = 1/2 in 
both 5D general relativity (obtained by setting d — A and a; — > 00 in (36), in agreement 
with Xu and Ma 2007) and in 4D Brans-Dicke with uj — 0. This is a refiection of the 
fact that a compact extra dimension can be represented as a scalar 'dilaton' field. A 
scalar field with a; = is exactly what is predicted in a such a situation (Carroll et al 
2002). In two or more extra dimensions the 00 parameter of the dilaton field becomes 
a function of the field and thus no longer behaves like a Brans-Dicke scalar (at least in 
the Jordan frame). 

The nonzero value of ^4 signals a breakdown of energy conservation in the four 
dimensional effective theory related to the potential $4, that is, related to pressure. 
We may understand this as arising from the fact that pressure is (by construction) 
isotropically distributed, even over dimensions which we cannot access. ^4 returns to 
zero when we take d = 3. A detection of a loss of pressure energy could then be used to 
unambiguously determine the true dimensionality of spacetime in Brans-Dicke gravity. 
However, given that pressure energy densities are so small compared to mass densities 
in any feasible experiment, there is little hope that we could measure (4 directly. 

It is also important to notice that f^, for example in (30), is a velocity norm over all 
dimensions, including the compact ones, whereas we can only see velocity components 
in the large dimensions. Thus motion in the compact dimensions will appear as an 
anomalous contribution to /iqo- This anomaly is eliminated if motion is restricted to 
the large dimensions. Such a restriction also returns ^4 to zero if pressure arises from 
microscopic particle velocities. 

On the other hand. Solar System experiments have succeeded in placing very strong 
constraints on the value of 7 which can in turn constrain the Brans-Dicke parameter u. 

Solar System experiments have set a limit a = I7 — 1| < 2.3 x 10~^ (Will 2006). 
Solving for uj gives 



We see that when d = 3, 00 ^ ±00 as a ^ 0. Given such tight bounds on a, we can 
constrain a; > 4 x 10^. Thus the scalar field is rather unimportant dynamically, if it 
exists at all. 

In contrast, when d > 3, cu is held very close to an order unity negative value of 
— {d—2)/{d—3). This confiicts with the original specification by Brans and Dicke (1961) 
of a positive cu. Assuming = for simplicity, we can easily derive the Hamiltonian of 
the theory. 



cu = 




(37) 




+ -V0 • V0, 




4a; (j) 
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where the conjugate momentum of the scalar field tt = 2c(;0/0. We see that if ui is 
negative, there is no ground state. In this way we confirm the well-studied result that 
the Jordan frame formulation of Brans-Dicke theory is unphysical (Faraoni et al 1999; 
Magnano and Sokolowski 1994; Cho 1992). Even for = 3, the presence of linear kinetic 
terms in the right hand side of (2) signals that the energy is not positive definite. The 
Jordan frame equations also fail to correctly specify the dynamics of the spin-2 graviton. 

These irregularities can be clarified by making the usual Einstein frame conformal 
transformation 

9t,u = Ht^u, (39) 



which, following a field redefinition 0/0o = exp[(0 — (po) / y uj + d{d — l)/4], casts the 
Lagrangian in the usual form of Einstein gravity with a minimally coupled free scalar 
field (p. This is the unique transformation which eliminates the problems mentioned 
in the previous paragraph (Magnano and Sokolowski 1994; Faraoni 1998). We achieve 
simplicity in the gravitational terms, however, at the expense of picking up a coupling 
between the scalar field and matter, thus violating the weak equivalence principle. 

If the Einstein frame is physical, experiments will be sensitive to the transformed 
metric Quiy. For our purposes, having already found the Jordan frame metric components 
to post-Newtonian order, we may simply multiply them by our expression (11) for the 
scalar field to effect the conformal transformation to appropriate order. Beginning with 
/loo to lowest order we find 

*o = 0o(l + j3Tf^)(-l + 2Ga) 

Since we know that Newtonian gravity holds to lowest order, we recognize a rescaled 
gravitational constant 

G = d-^/^+(f-f-G (41) 

d-l + {d-2)u ^ ^ 

so that hoo — 2GU . By a similar procedure, we find 

hi. = 2GU5i.- y^^^ — . (42) 

'^d-3/2 + (d-2)a; ^ ' 

The last factor in this expression is then recognized as our transformed PPN coefficient 

7- 

When (i = 3, 7 is identically 1, and we can get no information on uj. However, 
defining a = |7 — 1| and solving for u) gives 

3-rf+ (3/2-rf)a 

- = d-3 + (d-2)a ■ ^"^^ 

Thus we find that Solar System experiments constrain cu to a value of —1 for any value 
of d > 3 to within ±1.2 x 10~^n~^ where n = d — 3 is the number of extra dimensions. 
It is interesting to note that this is the same behavior predicted by string theory. The 
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gravitational terms of string theory effective actions take tlie form of Brans-Dicke gravity 
with uj = —1. Even after dimensional reduction, the predicted parameter value remains 
at —1 regardless of how many extra dimensions are considered (Lidsey et al 2000), 
analogous to our finding that c<j ^ — 1 for any value of > 3. It is also easy to derive 
the Einstein frame Hamiltonian and confirm that the theory with a; = — 1 is classically 
stable. 

It has been argued that if we consider the Einstein frame to be physical and if 
the minimal coupling of matter is taken as a fundamental principle, then we ought to 
effect the conformal transformation on the Lagrangian before adding in a matter term 
(Magnano and Sokolowski 1994). In this situation, however, the parameter uj has no 
independent influence outside of the transformed field 0, and we can learn nothing about 
it through the techniques discussed in this paper. 
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